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Abstract 

Background: Color traits in animals play crucial roles in thermoregulation, photoprotection, camouflage, and visual 
communication, and are amenable to objective quantification and modeling. However, the extensive variation in 
non-melanic pigments and structural colors in squamate reptiles has been largely disregarded. Here, we used an 
integrated approach to investigate the morphological basis and physical mechanisms generating variation in color 
traits in tropical day geckos of the genus Phelsuma. 

Results: Combining histology, optics, mass spectrometry, and UV and Raman spectroscopy, we found that the 
extensive variation in color patterns within and among Phelsuma species is generated by complex interactions 
between, on the one hand, chromatophores containing yellow/red pteridine pigments and, on the other hand, 
iridophores producing structural color by constructive interference of light with guanine nanocrystals. More 
specifically, we show that 1) the hue of the vivid dorsolateral skin is modulated both by variation in geometry of 
structural, highly ordered narrowband reflectors, and by the presence of yellow pigments, and 2) that the 
reflectivity of the white belly and of dorsolateral pigmentary red marks, is increased by underlying structural 
disorganized broadband reflectors. Most importantly, these interactions require precise colocalization of yellow and 
red chromatophores with different types of iridophores, characterized by ordered and disordered nanocrystals, 
respectively. We validated these results through numerical simulations combining pigmentary components with a 
multilayer interferential optical model. Finally, we show that melanophores form dark lateral patterns but do not 
significantly contribute to variation in blue/green or red coloration, and that changes in the pH or redox state of 
pigments provide yet another source of color variation in squamates. 

Conclusions: Precisely colocalized interacting pigmentary and structural elements generate extensive variation in 
lizard color patterns. Our results indicate the need to identify the developmental mechanisms responsible for the 
control of the size, shape, and orientation of nanocrystals, and the superposition of specific chromatophore types. 
This study opens up new perspectives on Phelsuma lizards as models in evolutionary developmental biology. 

Keywords: Physics of biology. Structural colors. Pigmentary colors. Color patterns, Iridophores, Chromatophores, 
Reptiles, Lizards, Phelsuma 
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Background 

Vertebrate skin coloration provides a promising model 
system for exploring the link between genotypes and phe- 
notypes in an ecological and phylogenetic framework. 
Indeed, color traits play crucial roles in thermoregula- 
tion, photoprotection, camouflage, and visual communi- 
cation [1-4], and can vary extensively among and within 
species and populations. Moreover, colors and color patterns 
are amenable to objective quantification and modeling, pro- 
viding an opportunity for integrated analyses of phenotypic 
variation. In particular, non-mammalian vertebrates, for 
example, squamates (lizards and snal<es), exhibit a broad 
range of pigmentary and structural colors, generated by 
different types of chromatophores. In addition to melano- 
phores, which produce black/brown melanins, squamates 
develop xanthophores and erythrophores, containing yellow 
and red pigments, respectively. These pigments are typically 
either pteridines, which are synthesized in situ from guano- 
sine triphosphate, or carotenoids, which are metabolized 
from food in the liver and transported to skin via the cir- 
culatory system [5,6]. Squamates additionally possess 
iridophores, which do not contain any pigment but gene- 
rate structural coloration through interference of light 
waves with transparent guanine nanocrystals [7-9]. The 
spatial arrangement of all these cell types generates a 
broad range of colors and color patterns. 

Despite such a high potential for complexity and diver- 
sity, lizards and snakes remain relatively under-represented 
in evolutionary developmental studies in general [10], and 
in analyses of color-pattern evolution in particular (for ex- 
ample, in comparison with other vertebrates or insects) 
[3,4,11,12]. Although chromatophores have been described 
in some squamates [8,9,13-16], and their melanin pathway 
has been associated with adaptive color variation [17,18], 
few data are available on the mechanisms that generate ex- 
tensive variation in non-melanic pigments and structural 
colors in this lineage. The presence of such a variety of 
colors in squamates makes them appropriate models for 
investigating the essentially unknown genetic, develop- 
mental, and physical mechanisms generating a diversity of 
phenotypes through interactions between different types 
of chromatophores. 

Here, we used an integrated approach that combines 
histology, Raman and UV spectroscopy, mass spectro- 
metry (MS), optics, and numerical simulations to investi- 
gate the morphological basis and the physical mechanisms 
generating variation in color traits in five representative 
tropical day gecko species of the genus Phelsuma, More 
specifically, we show that 1) the reflectivity of the white 
belly, and of the dorsal pigmentary red spots and 
stripes, is increased by underlying, structurally disorga- 
nized broadband reflectors, whereas 2) the hue of the 
vivid dorsolateral background coloration is generated 
by structural, highly ordered narrowband reflectors. 



and is modulated both by the photonic crystal geometry 
and by a layer of yellow pigments. Most importantly, we 
show that these interactions require precise colocalization 
of red and yellow pigment cells with iridophores charac- 
terized by different organizations of nanocrystals. 

Results and discussion 

Tropical day geckos of the genus Phelsuma provide excel- 
lent models to study the morphological basis of color- 
pattern variation. Most of the approximately 45 currently 
recognized species [19,20] are characterized by an off-white 
ventral and a vivid dorsolateral coloration, with a blue/ 
green background and red (sometimes brownish) spots or 
stripes on the back (Figure la). This vibrant and contrasting 
color scheme is likely to be involved in mate recognition 
and/or camouflage [1,2], and is highly variable within 
and among species. Not only can the background color 
differ substantially among individuals (as, for example, in 
Phelsuma grandis; Figure lb, bottom panels), but also 
the size, shape and intensity of the red dorsal markings. 
In addition, several species, including Phelsuma klemmeri, 
Phelsuma quadriocellata and Phelsuma lineata, have la- 
teral black/brown spots or stripes, and a few, such as 
P, klemmeri lack red markings and exhibit a blue/brown 
background (Figure lb, top panel). 

The skin organization in Phelsuma geckos (Figure 2a) 
resembles that of other lizards [9]: chromatophores are 
absent from the thin epidermis, but are abundant in the 
thick dermal layer, which contains, from top to bottom, 
yeUow xanthophores or red erythrophores (present in 
green and red skin respectively), iridophores containing 
nanocrystals, and dark-brown melanophores. Here, we 
report on the chemical and optical analyses of pigments 
and guanine crystals, and provide a multflayer optical 
model that describes the interactions among different 
types of pigmentary and structural chromatophores ge- 
nerating a variety of skin colors in Phelsuma geckos. 

Melanophores generate dark lateral spots and stripes 

The skin of all geckos studied here, regardless of color, 
contains two types of melanophores (Figure 2a): smaU 
light-brown ceUs that, when associated with iridophores, 
form 'dermal chromatophore units' [21], and large dark- 
brown cells at the base of the dermis. The former are 
dendritic ceUs with processes extending through the 
layer of iridophores for translocation of melanin gra- 
nules in response to hormones, resulting in darkening of 
the skin [9,22]. For example, the melanophore processes 
cover the iridophores in the dorsal skin of P. klemmeri, 
giving it a light-brown appearance. Moreover, in dark-brown 
lateral spots and stripes of P, klemmeri, P, quadriocellata, 
and P, lineata, small melanophores are abundant on top 
of iridophores, that is, in the layer typically occupied by 
xanthophores or erythrophores in other parts of the body 



Saenko et al. BMC Biology 201 3, 1 1 :1 05 
http://www.bionnedcentral.conn/1 741 -7007/1 1 /1 05 



Page 3 of 1 2 



(a) 




Species 



Background 



(b) 



P. kfemmeri 



Background 



Pattern 



— E 



— c 







dorsal 


lateral 








• » 1 


m 


— 






• — 






P. standingf 




— 1 






Pr serrattcsuda 




• — 






P. kl^mmerl 










p. antanosy 




• — 






P, quadrioceHata 






• 




p. q. cf, bimsculata 




• — 


• 




Pr parva 


1 — II — II — 1 ■ H 




9' 




P. pii$iit0 










P. keiy 




• 






P, tineata 










p. comarensis 




• — 




— 


P. dorsMttata 






• 














P, Istfcauda 










p. pssteuri 










P. robsitmentensl 










Pr v-nigra 




• 
















P. hieischerf 










P. bergtiofi 


□ ■□BB 








P. matamakitKf 


UUUHLJ 










□ □□■M 






: = 


P. aubia 










P. rsvenata 










P. nigristn&ts 










P. modesta 








— 












P. barboufi 




• — □ 


O 




if. fironki 




... — □ 
















P. gijentherf 










p. ros^guf^fius 










f^j, iJkJf ilLffrti-a 


1 — 1 H L — 1 H HI 


• — ■ 






P. cepediana 










P. inexpectata 






























p. brevicepS 










P. borai 










mi^abflis 




• — □ 


o 














P. vanheygeni 










:P. astiiata 




















P. guttata 










P. madagascariensis 










P. seippi 










P. kochi 










P. grandis 




• — 1 






p. parkeri 










P, abbot! 




• — 








P. Hneata 



Dorsal pattern 




P. grandis 




Ventral 



Figure 1 Color and pattern variation in tlie genus Phelsuma. (a) Currently recognized species [19,20] witli corresponding ranges of bacl<ground 
coloration (indicated by filled squares), varying from gray/brown to yellow/green to turquoise-blue, and dorsal and lateral patterns (indicated for 
illustration on Phelsumo lineoto with an arrow and an arrowhead, respectively). Empty shapes indicate absence of the corresponding color or pattern. 
Dorsal patterns vary in both color (typically different shades of red) and shape (spots, stripes, and transverse bars). Lateral patterns appear as spots 
(sometimes surrounded by a blue ring) or broad stripes. Species indicated in bold were used in this study, (b) Skin samples represent variation found 
across the genus. The skin on the belly (shown here for Phelsuma grandis, lowest panel) is off-white in the majority of species. 
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Figure 2 Pigmentary and structural colorations in Phelsuma geckos, (a) Semi-thin cross-sections of sl<ins of different colors. Two types of 
melanopliores are indicated (arrows and arrowlieads, respectively), and iridophores are outlined with dashed lines. Ep, epidermis; X, xanthophores; 
Er, erythrophores. Bar = 10 pm. (b) Pteridin pigments were removed with NH4OH (here in Phelsumo grandis, individual number 3), revealing the 
remaining structural color produced by the iridophores. Bar = 0.2 mm. (c) Red pigments in dorsal markings of Phelsumo quodriocelloto and 
Phelsumo lineoto can change color when the pH of the Ringer solution is lowered or when an oxidant (NaN02) is added, respectively. Bar = 0.2 
mm. (d) Representative electron micrographs of iridophores in skin of different colors. Bar = 1 pm. Note the highly disordered guanine crystals in 
the white and red skin, (e) Mechanical pressure and dehydration (here applied to the green skin of P. grondis, individual number 2 after removal of 
the yellow pigment) lead to a blue shift of structural green (for supplementary movies, see Additional file 3; see Additional file 4). 
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(Figure 2a). In P. lineata, melanophores are found in 
combination with red erythrophores in the red-brown re- 
gions of the lateral stripes. Hence, in addition to their 
likely involvement in darkening of the skin, the small me- 
lanophores form dark lateral spots/stripes in some 
Phelsuma geckos, whereas the larger and darker melano- 
phores might protect the deeper layers from harmful UV 
radiation [22]. However, the homogeneous distribution of 
these two types of chromatophores in the skin of any 
color suggests that neither of them substantially contribute 
to color variation in Phelsuma geckos. 

Pteridins contribute to red and green colors of 
dorsolateral skin 

We could dissolve in NH4OH (but not in acetone) both the 
red and yellow pigments, found respectively in red-skin 
erythrophores and green-skin xanthophores (Figure 2b). 
This indicates that both pigments belong to the pteridine 
class rather than to the carotenoid class [23,24]. We 
confirmed these findings by Raman spectroscopy and 
UV spectroscopy/MS (see Additional file 1). The Raman 
spectra (see Additional file 1: Figure SI) of both yellow 
xanthophores and red erythrophores in P, quadriocellata, 
P. lineata, and P. laticauda are similar to that of 
xanthopterin [25], a known pigment in squamates [15]. 
Conversely, although the Raman spectrum of yellow 
xanthophores in P. grandis also indicates xanthopterin, the 
spectrum obtained for the red erythrophores in this spe- 
cies is significantly different, suggesting the presence of 
other molecules. 

We then performed UV chromatography/MS on pigments 
extracted from the red and green sldn of P, quadriocellata, 
P. lineata, and P. grandis, and compared their spectra 
with those of available pteridin standards. We confirmed 
the presence of xanthopterin in the yellow chromato- 
phores of P. grandis, and in both the yellow and red 
chromatophores of P. lineata. We also identified biopterin 
in the yellow and red chromatophores of P. quadriocellata 
and in the red chromatophores of P, grandis, and identi- 
fied sepiapterin in the yellow chromatophores of P, 
grandis. Finally, we found three unidentified molecules 
(probably pterins, based on their UV absorption spectra 
(see Additional file 1: Figure S2)) in all three species. It is, 
however, unclear which components contribute most sub- 
stantially to the final pigmentary colors of the skin. 

Remarkably, although pigment compositions differ sub- 
stantially between P, lineata and P, quadriocellata, the pig- 
ment compositions of xanthophores and erythrophores 
are identical within a species (see Additional file 1: Figure S2a), 
strongly suggesting that the actual in vivo colors of these 
chromatophores (red or yellow) are determined by some 
other factors, such as the pH of the cellular environment of 
the pigment (as is the case for plant anthocyanins [26]) or 
the redox state of the pigment molecules (as has been 



shown for ommochromes in dragonflies [27]). Supporting 
this hypothesis, the red color of dorsal sldn changed to yel- 
low when we lowered the pH of, or added an oxidant 
(NaN02) to, the skin of P, quadriocellata or P, lineata, re- 
spectively (Figure 2c). Even though the exact in vivo pro- 
cesses generating variation in the pH or redox states are 
unknown, these results indicate that these mechanisms pro- 
vide an additional source of color variation in squamates. 

Ordered multilayer interference reflectors generate 
structural blue or green 

In all regions of the body, irrespective of its color, the skin 
of Phelsuma geckos contains iridophores, with guanine 
nanocrystals, either well- organized in parallel layers (in 
the light-brown dorsal skin of P. klemmeri and in the 
green/blue skin of all other species), or highly disorga- 
nized (in the white and red skin of all animals examined 
here) (Figure 2d). Using transmission electron microscopy 
(TEM) images, we measured the orientation of crystals 
relative to the skin surface, and computed the full-width 
half maximum (FWHM) of the Gaussian distribution of 
crystal orientation, and the ratio between the amplitude A 
of the Gaussian curve and the amplitude jo of the back- 
ground generated by randomly oriented crystals (see 
Methods; see Additional file 2: Figure S3d). The A/jq ratio 
represents the relative amount of crystals oriented parallel 
to the skin surface and varies from zero (when A = 0 in a 
purely random system) to infinity (when A is high and jo 
is close to zero in a perfectly ordered system). Table 1 
shows that blue/green skin is characterized by a narrower 
angular distribution of nanocrystals and a higher A/jq ratio 
than white or red skin. This indicates that the ordered 
nanocrystals generate an intense narrowband structural 
blue or green color through a multilayer interference 
mechanism. This result is consistent with the pigment- 
removal experiments (Figure 2b) discussed above, revealing 
the structural blue color of iridophores that are typically 
covered by the yellow pigmentary color of xanthophores in 
the green skin of some P, grandis individuals. Note that 
crystal spacing (hence, structural color) varies both among 
and within individuals (Table 1). 

To further test experimentally if the multilayer is respon- 
sible for the structural color component of blue/green skin 
in Phelsuma, we applied mechanical or osmotic pressure to 
modify the distance between the guanine nanocrystal layers 
and therefore induce changes in structural coloration [16]. 
We found that mechanical pressure or dehydration applied 
to de-pigmented green skin (hence, exhibiting structural 
green) resulted in a shift to structural dark blue (Figure 2e; 
see Additional files 3 and 4: supplementary movies). Note 
that complete dehydration made the layer of iridophores 
transparent, revealing the underlying melanophores 
(Figure 2e), because a severe decrease in crystal spacing 
either obliterates coherent interferences (that is, the 
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Table 1 Geometric parameters of guanine crystals in Phelsuma species 

Species Green/blue White Red 



Orientation Crystal Crystal Orientation Crystal Crystal Orientation Crystal Crystal 

A/yo (FWHM)^ height, nm'' spacing, nm'' A/yo (FWHM) height, nm length, nm A/yo (FWHM) height, nm length, nm 



P. grand is (1) 


15.0 (32°) 


80.8 ± 


12.2/ 


96.0 ± 20.8/ 


2.7 


(72°) 


84: 


h44 


173 d 


: 84 


1 .3 (82°) 


68 ± 


: 34 


128 


± 70 






78.5 ± 


10.9' 


104.7 ± 25.5' 
























P. grondis (2) 


9.7 (49°) 


81.0 ± 


12.9 


84.9 ± 22.7 


2.0 


(61°) 


78 : 


h 39 


169 d 


: 96 


1.1 (53°) 


67 ± 


: 33 


124 


± 82 


P. grondis (3) 


13.0 (47°) 


79.4 ± 


10.7 


99.2 ± 10.6 


2.6 


(59°) 


75 : 


h 36 


171 d 


: 91 


1 .7 (46°) 


71 ± 


: 37 


130 


± 73 


P. 1 at i Cauda 


13.3 (39°) 


82.5 ± 


11.4 


97.3 ± 12.6 


2.2 


(60°) 


71 : 


h41 


194 d 


: 140 


2.6 (61°) 


75 ± 


: 42 


154 


± 87 


P. quadriocellata 


20.3 (40°) 


70.0 ± 


22.7 


37.9 ± 13.8 


6.4 


(49°) 


84: 


h42 


209 d 


: 126 


4.8 (67°) 


69 ± 


: 41 


171 


± 103 


P. liiemmeri 


9.6 (32°) 


68.0 ± 


10.8 


93.2 ± 26.3 


4.5 


(58°) 


72 : 


h47 


265 d 


: 209 


NA^ 


NA^ 




NA^ 




P. lineata 


6.9 (42°) 


74.0 ± 


11.1 


109.6 ± 21.8 


4.8 


(50°) 


94 : 


h 57 


226 d 


: 137 


2.1 (58°) 


68 ± 


:43 


170 


± 107 



^The ratio A/yo and the full-width half maximum (FWHM) of the Gaussian curve describe the crystal orientation relative to the skin surface (see Methods). 

^Mean ± standard deviation are given for crystal height and spacing between layers of well-organized crystals in green and blue skin, and for crystal height and 

length in white and red skin. Note that the parameters describing crystal orientation differ substantially between blue/green and white/red skin. 

"^The two numbers shown for crystal height and spacing of green/blue skin in P. grandis individual 1 correspond to samples taken from the back and the neck 

(blue and yellowish-green, respectively; see Figure 4b). Note that the spacing between crystal layers is substantially larger in the neck sample. Additional data are 

provided for blue and green skin samples, including a comparison of the simulated colors to the real colors seen in the skin after pigment removal 

(see Additional file 2: Table S3). 

'^NA, data not available because of the absence of red skin in P. klemmeri. 



entire stack of crystals behaves as a single block of gua- 
nine), or causes a shift in the reflectivity peak beyond 
the optical visible range, that is, below 390 nm. 

Optical simulations of ordered multilayer Interference 
match the reflectivities and colors produced 
by iridophores 

We used a multilayer model (Methods; see Additional 
file 2) to simulate the reflectivity and colors produced 
by organized crystals using parameters such as crystal 
height (that is, thickness) and spacing between crystal 
layers, estimated from TEM images collected for a se- 
lection of green and blue skin samples (Table 1; see 
Additional file 2: Table S3). The measured and modeled 
reflectivities of four skin samples, covering the range of 
structural colors seen in the Phelsuma individuals stu- 
died here, are shown (Figure 3a). Although we did not 
have TEM data available for the sample of P. grandis 
number 4, its color was similar to that of the P, grandis 
number 1 neck sample. In addition, we compared the 
RGB (red, green, blue) colors predicted by the model to 
color pictures of all samples after pigment removal (see 
Additional file 2: Table S3). Overall, the reflectivities and 
the RGB colors predicted by our model were similar, but 
not always identical, to those perceived on the animals. 
This small discrepancy can be explained by, for instance, 
changes in cytoplasm osmolarity (and hence, the spacing 
between crystal layers) during skin preparation for TEM 
[14-16]. In addition, we found high reflectivity in the UV 
spectrum (Figure 3a), which might be produced by a series 
of mechanisms such as coherent scattering on arrays of 
collagen fibers in the dermis [28] or incoherent scattering 
[29] on, for example, melanosomes or guanine crystals. 



The collaborative interplay between the structural and 
pigment components is illustrated (Figure 3b); we mea- 
sured the reflectivity of green skin before and after yellow 
pigment removal (green and blue solid lines, respectively). 
The modeled reflectivity for an arbitrarily chosen crystal 
size (70 nm) and spacing (mean = 30 nm, SD = 13 nm), 
shown as dashed blue line, closely resembles the measured 
reflectivity of green skin after pigment removal. When a 
3 (im yellow pigment layer is included in the model, the 
resultant reflectivity (Figure 3b, dashed green line) is a 
good match to the measured reflectivity of green skin be- 
fore pigment removal. A similar reflectivity distribution 
(red line) was obtained as the direct product of the blue 
line with the measured normalized transmittance of the 
yellow pigment (orange dashed line). Hence, the green 
color of the skin is due to the structure-based blue reflec- 
tion filtered by the yellow-transmitting pigment. 

Optical systems based on multilayer interference usu- 
ally show a spectral shift in reflectivity as the angle of the 
incident light changes [30]. This angular dependence is 
determined by the refractive indices of the alternating 
layers, and the ratio between the optical paths at different 
angles (that is, path length x the refractive index of the 
medium) [31,32]. Under very particular combinations of 
the refractive indices and thicknesses of its two compo- 
nents, a multilayer reflector can exhibit greatly reduced 
angular dependence [33]. Only one of these conditions is 
met in Phelsuma geckos: the optical paths in low {ncd^ 
and high {n^^ index layers are close, making the multi- 
layer a narrowband optimal reflector as a result of the 
collaborative effects of crystal and spacing interferences. 
However, the iridescence in Phelsuma skin is much 
weaker than that predicted by the model (see Additional 
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Figure 3 Experimental and modeled reflectivities of Phelsuma skin, (a) Measured skin reflectivities (solid lines) after removal of pigments 
(skin colors varied from deep-blue to yellowish-green) compared with modeled reflectivities (dashed lines) based on crystal size and spacing 
(Table 1). For Phelsuma grandis number 4, the crystal geometry parameters were taken from the P. grandis number 1 neck sample that exhibited 
a similar structural color. Note that the UV peaks in the measured reflectivities are probably caused by scattering on melanosomes, on iridophore 
crystals, or on dermal collagen fibers, (b) Normalized reflectivity of green skin before and after yellow pigment removal (green and blue solid 
lines, respectively). Modeled multilayer responses for a crystal size of 70 nm and a spacing of 30 nm (assuming a standard deviation of 13 nm) 
are also shown with (dashed green line) and without (dashed blue line) a 3 pm thick pigment layer on top. The direct product of structural blue 
reflectivity with normalized yellow pigment transmittance (orange dotted line) generates the plain red line, confirming the mechanism of 
structural color filtering by the top pigment layer, (c) Reflectivity measured on P. grandis white skin (black line) and on red skin before (red solid 
line) and after (red dashed line) red pigment removal. Reflectivity intensities of ordered and disordered iridophores are comparable. 



file 5: Figure S5a) because of another reason, namely, the 
variable orientation of iridophores in the skin reduces or 
annihilates iridescence by averaging the contributions 
from different cells (see Additional file 5: Figure S5b). 

Incoherent scattering by disorganized crystals enhances 
the reflectivity of red dorsal and white ventral skin 

Guanine crystals in iridophores of green and blue skin 
are well-organized, and therefore contribute extensively 
to the final background color of Phelsuma lizards. Re- 
markably, iridophores are also abundant in both the white 
ventral skin and red dorsal markings, but these cells are 
characterized by crystals with mostly random orientations 



and broad size distribution (Figure 2d; Table 1). This gives 
rise, as in some insects and fish [34-39], to incoherent 
scattering. In other words, in contrast to iridophores in 
blue/green skin, which generate a narrowband distribution 
of reflected frequencies, iridophores in white and red skin 
form broadband reflectors with overall reflectivity of simi- 
lar intensity to that measured on coherently scattering 
structures (Figure 3c). Hence, these cells with disorganized 
guanine crystals not only produce the white color of the 
belly, but also significantly enhance the reflectivity, and 
therefore the visibility, of the red pigmentary patterns on 
the back. In the absence of iridophores, the red spots would 
appear less bright because much of the incident light would 
be absorbed by the underlying tissues. The mechanisms 



Saenko et al. BMC Biology 201 3, 1 1 :1 05 
http://www.bionnedcentral.conn/1 741 -7007/1 1 /1 05 



Page 8 of 1 2 



responsible for the distribution of ordered versus disordered 
crystals and for this spectacular colocalization of pigment 
cells with a particular type of iridophores are unknown. 

Variation in skin coloration is explained by both 
structural and pigmentary components 

Variation in background coloration and dorsolateral pat- 
terns in Phelsuma geckos (Figure 1) can be generated by 
changes in pigmentary and structural components. For 
example, variation in the intensity of the red dorsal pat- 
terns can be caused by differences in the density of 
erythrophores in the top dermal layer and/or in concen- 
trations of red pigment in these cells. The color variation 
obtained from a simulated red pigment layer of variable 
thickness (0.1 to 4 (im) on a fully reflecting background 
was a good match to the variation in intensity of the red 
markings seen among species or individuals (Figure 4a). 
Similarly, different densities of xanthophores and/or con- 
centrations of yellow pigment, in combination with struc- 
tural colors produced by iridophores, can lead to variation 
from blue to green. For example, the background skin 
color of P. grandis individual number 1 (Figure 4b) varies 
from a greenish-blue to yellowish-green depending on the 
area of the body involved. This variation is explained by 
an increase of about 10% in spacing between guanine crys- 
tals (96.0 nm in the back skin versus 104.7 nm in the neck 
skin; Table 1). The colors of these two different zones on 
the skin (Figure 4b) are in agreement with the predictions 



of our optical model (Figure 4c). Note that similar vari- 
ation in skin color could also be due to modification of 
the thickness or density of the yellow pigment layer with- 
out any change in the structural component (Figure 4d). 

Conclusions 

Extensive variation in background coloration, and in 
dorsal and lateral color patterns, is present within and 
among species of the genus Phelsuma (Figure 1). Our 
analyses indicate that this variation is generated by a 
combination of features associated with pigmentary and 
structural color chromatophores (as previously reported 
in birds [40,41]). The black lateral spots and stripes 
found in some species, and the light-brown background 
coloration of P, klemmeri, are due to melanophores oc- 
cupying the upper layers of the dermis. Blue skin color 
is solely due to iridophores with well-organized guanine 
nanocrystals (that is, narrowband multilayer interference 
reflectors), whereas green skin is produced, depending 
on the species and the individual, either by structural 
green or by the interaction of structural blue with yellow 
pigments (xanthophores). Dorsal marks are formed by 
red erythrophores, the reflectivity of which is enhanced by 
iridophores with disorganized crystals (that is, broadband 
reflectors producing incoherent scattering). In addition, 
the hue produced by erythrophore and xanthophore pig- 
ments is pH-dependent or redox-dependent in some spe- 
cies. Most importantly, we show that the color patterns of 
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Figure 4 Colors simulated with the multilayer model, (a) Red colors simulated with varying thickness (0.1 to 4 pm) of a red pigment layer on 
top of a white reflector and comparison to red markings of different animals, (b) Phelsuma grondis individual number 1. (c) Simulated colors produced 
by a 7 |jm yellow pigment layer on top of a multilayer interference reflector with varying spacing (40 to 105 nm) between layers of crystals 80 nm 
thick. Double arrows indicate the spacing measured between crystal layers on the dorsal and neck skin of the individual, (d) Simulated colors 
produced with varying thicknesses (0 to 7 |jm) of a yellow pigment layer on top of a blue reflector. 
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Phelsuma always require precise colocalization of different 
sets of interacting pigmentary and structural elements. 
For example, yellow and red chromatophores are associ- 
ated with iridophores with ordered and disordered 
nanocrystals, respectively. This indicates the need to iden- 
tify the developmental mechanisms responsible for the 
superposition of specific chromatophore types, opening 
up new perspectives for Phelsuma lizards as model orga- 
nisms in evolutionary developmental biology. 

We also show that the exact combination of parame- 
ters producing a given color is difficult to predict with- 
out experimental estimation of parameter values (such 
as, through electron-microscopy and spectroscopy ana- 
lyses). Hence, exploring the genetic and developmental 
bases of phenotypic variation in ecologically important 
color traits in squamate reptiles will require integration 
of the structural and pigmentary color components and 
their interactions. The fact that red/yellow coloration in 
Phelsuma is based on pteridines, the synthesis of which 
is controlled by enzymes that have been well studied in 
model organisms [42], readily suggests candidate genes 
for this pigmentary aspect of color variation. Further- 
more, although guanine crystal formation in iridophores 
is poorly understood, and the genetic determinants of 
variation in size, shape, and orientation of nanocrystals 
are virtually unknown, several enzymes associated with 
aberrant iridophore phenotypes [43,44] are candidates 
for variation in structural coloration in Phelsuma, 

Methods 

Maintenance of, and experiments on, animals were 
approved by the Geneva Canton ethical regulation autho- 
rity (authorization 1008/3421 /IR) and performed in ac- 
cordance with Swiss law. 



Animals and phylogenetic mapping 

Adult P, klemmeri (n = 1), P, quadriocellata (n = 1), P, lineata 
(n = 2), P, laticauda (n = 1), and P, grandis (n = 5) were 
obtained from the pet market or bred in our laboratory. 
On the phylogenetic tree of the genus Phelsuma [19,20], 
we mapped (Figure la) the presence or absence of the fol- 
lowing phenotypic characters (based on previous reports 
and our own observations [20,45-50]): 1) dorsolateral back- 
ground coloration; 2) red/brown dorsal pattern; and 3) 
black/brown lateral pattern. Most Phelsuma geckos exhibit 
a vivid dorsolateral coloration, with a background typically 
of one or more colors (for example, bright green, yellow, 
or turquoise-blue, but sometimes dull gray or brown). Al- 
most all species have distinct dorsal marks consisting of 
red to brown spots and stripes of various shapes. In 
addition, some species have black/brown lateral stripes or 
spots. 



Skin histology and TEM 

Skin samples were placed in Ringers solution for micros- 
copy observation. Cross-sections of 14 to 16 (im were pre- 
pared from skin embedded in optimum cutting temperature 
compound on a Leica (Wetzlar, Germany) CM1850 cryo- 
stat. For TEM, skin pieces of approximately 1 cm^ in size 
were frjced overnight at 4°C in 2% glutaraldehyde and 4% 
paraformaldehyde in 0.05 mol/L sodium cacodylate buffer 
(pH 7.4), rinsed with 0.1 mol/L cold sodium cacodylate, 
and frjced for 1.5 hours on ice in the same buffer with 
the addition of 1% osmium tetraoxide. Samples were then 
rinsed in water, stained with 1% uranyl acetate for 2 
hours, gradually dehydrated in ethanol, rinsed three times 
in propylene oxide, incubated overnight in 1:1 propylene 
oxide/resin, and finally embedded in epoxy resin (Epon). 
Semi-thin (1 to 2 (im) and ultra-thin (80 to 90 nm) 
cross-sections were cut with a diamond knife on a Leica 
UCT microtome. Semi-thin sections were examined under a 
light microscope. Ultra-thin sections were placed on 
formavar-coated grids, and viewed with a Tecnai™ G2 Sphera 
(FEI) TEM at 120 kV to visualize intact guanine crystals, be- 
cause post-staining inevitably results in loss of crystals [13]. 
The grids were then post-stained with uranyl acetate and 
lead citrate, and viewed with the TEM again (in this case, 
spaces that were occupied by crystals appear white). 

Optical modeling of skin reflectivity and color 

For each individual, 10 to 20 TEM images of stained sec- 
tions of iridophores of white, red, and green skin were col- 
lected at low (x 1,500) magnification, and 20 to 30 images 
of unstained sections were collected at high (x 25,000) mag- 
nification for green/blue skin. To analyze the size and 
orientation of crystals relative to skin surface, white 'holes' 
(corresponding to guanine crystals dissolved by uranyl acet- 
ate/lead citrate during post-staining) on low magnification 
images were fitted with ellipses (in JMicro Vision [51]), 
whose geometric parameters (length, height, and orienta- 
tion of major axis) were computed subsequently. 

The Gaussian distribution of crystal orientation was 
calculated as: 

where is the intersection of the curve with the y-axis 
(that is, the background level of randomly oriented crys- 
tals); A is the amplitude of the Gaussian curve; 6 is the 
angle of the smallest axis of the ellipse relative to the nor- 
mal to the skin surface; Oc is the central value of the distri- 
bution, and is equal to 0; and iv, the standard deviation, is 
related to the FWHM of the Gaussian curve, defined as 

2v^2 X ln(2) X w. 

The A/yo ratio varies from zero (when A = 0 in a 
purely random system) to infinity (when A is high and 
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yo is close to zero in a perfectly ordered system) as sche- 
matically illustrated (see Additional file 2: Figure S3d). 

For optical modeling in blue/green skin, the thickness 
of crystals and the spacing between well-aligned crystal 
layers were measured on images of unstained sections 
taken at x25,000 magnification (Table 1; see Additional 
file 2: Table S3). These values were used in a multilayer 
model (see Additional file 2) that simulates the reflecti- 
vity and color produced by alternating layers of guanine 
crystals and cytoplasm, each with a different thickness 
and refractive index [14,52,53]. 

To account for the contribution of pigments to skin 
reflectivity and hue, a top layer of pigment with variable 
thickness and wavelength-dependent absorption was in- 
cluded in the model. The optical transmission of the pig- 
ments measured on skin cryosections was fitted to a 
Drude-Lorentz model and included in the multilayer 
(see Additional file 2). As the intensity of light reflected 
coherently by organized photonic structures is much 
more intense (reflectivity can be close to 1) than that 
reflected by incoherent scattering, the latter was not 
taken into account. 

Pigment chemical analysis 

To determine pigment composition in xanthophores and 
erythrophores, skin samples were treated with 100% 
acetone and 30% NH4OH, known to dissolve carotenoids 
and pteridins, respectively [23,24]. MS was performed on 
pteridin standards (neopterin, isoxanthopterin, xanthopterin, 
pterin, 6-biopterin and sepiapterin) and on extracted 
pigments by UHPLC/DAD/ESI-QTOF (Agilent Tech- 
nologies 1290 Infinity, ESI-QTOF model G6530A). A 
Poroshell C8 analytical column (100 x 3.0 mm internal 
diameter; particle size 2.7 [im) was used for the sepa- 
ration. The column oven temperature was set at 40°C. 
The binary mobile phase consisted of 0.4% acetic acid 
in water (solvent A) and acetonitrile (solvent B), and 
the flow rate was set to 0.5 ml/min. UV/visible spectra 
were recorded between 190 to 600 nm, and the injection 
volume was 1 [A, Mass spectra were acquired in positive 
and negative mode using electrospray ionization with ni- 
trogen as the nebulizing gas, and recorded for the range of 
m/z 50 to 1700. Drying gas flow was 10 1/min, with a frag- 
mentation voltage of 120 V, drying gas temperature of 
300°C, nebulizer pressure of 40 psi, and capillary voltage 
of 3500 V. Compounds were identified with MassHunter 
Qualitative Analysis Software (Agilent Technologies) by 
analysis of their UV and high-resolution MS spectra (see 
Additional file 1). 

Raman spectroscopy was performed on cryosections of 
red and green/yellow skin. This non-invasive technique 
[54] determines the chemical structure of a target mo- 
lecule by measuring the energy shift of an incident laser 
light after its interaction with the molecule. The bespoke 



micro-Raman system was composed of a half-meter focal 
length spectrometer coupled to a nitrogen-cooled Prince- 
ton CCD detector. The excitation source was an argon 
laser with a wavelength of 514.5 nm. The full collection 
of Raman spectra on different pigments is shown (see 
Additional flle 1). 

Optical measurements 

For extended UV reflectivity measurements (300 to 
800 nm), a probe (QR400-7-VIS/BX; Ocean Optics) was 
connected to the scanning Xenon source of a spectro- 
scopic Woollam ellipsometer with silicon diode as de- 
tector. Lock-in detection was used for ambient light 
noise reduction. For monitoring color change in real 
time (see Additional files 3 and 4, movies) and for im- 
aging of the skin samples (see Additional file 2: Table S3), a 
Thorlab color CCD camera was mounted on a Leica Ml 25 
stereo microscope with a high field depth. With CCD, a 
systematic underestimate of RGB numbers of about 5% 
(that is, darker colors) was seen, but the color balance was 
rendered correctly. During the mechanical pressure and 
dehydration experiments, skin color was measured in situ, 
and compared with simulations in which the distance 
between the crystals was the only adjusted parameter. 

Additional files 



Additional file 1: Pigment analyses. Results of Raman and UV/mass 
spectroscopy on pigments. 

Additional file 2: Optical model. Description of the optical model used 
to simulate skin colors. 

Additional file 3: Movie of mechanical pressure in Phelsuma 
species. Real-time movie showing how pressure manually applied with 
tweezers locally modified the structural color, and how the relaxation 
process drove the skin back to its original color. Note that the yellow 
pigments had been removed before the experiment. 

Additional file 4: Movie of osmotic pressure in Phelsuma species. 

Accelerated movie (xl6) showing the effect of dehydration (in air) and 
rehydration on the structural color generated by iridophores. Note that 
the yellow pigments had been removed before the experiment. 

Additional file 5: Angular dependence. Angular dependence of the 
multilayer reflectivity in green skin of Phelsuma grandis. 
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